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1. 


INTRODUCTION 


The  Naval  Environmental  Prediction  Research  Facility  (NEPRF)  has 
collected  a data  base  of  upper-air  soundings  recorded  by  radar  picket 
ships  off  the  west  coast  of  the  United  States;  the  data  base  relates  to  the 
period  October  1956  through  June  1965,  and  the  area  29°-51°N,  120°-140°W. 
Under  the  terms  of  Contract  No.  N00228-76-C-3139,  this  data  base  is  to 
be  used  to  develop  a synoptic  climatology  of  refractive  conditions 
significantly  affecting  the  tropospheric  propagation  of  electromagnetic 
waves . The  development  is  to  be  carried  out  in  accordance  with  the 
following  three  Tasks: 

Task  One:  Monthly  climatology  of  refractive  conditions 

a.  Define  for  each  month  the  mean  number  of  days  during 
which  given  refractive  structures , thicknesses  and 
intensities  can  be  expected. 

b.  Define  a monthly  measure  of  variability  (variance  and 
standard  deviations)  for  each  of  the  indicated  conditions. 

Task  Two:  Persistence  of  refractive  conditions,  and  transition 

probabilities 

a.  Determine  mean  number  of  days  for  which  a given 
refractive  condition  can  be  expected  to  persist  and  the 
variability  of  such  persistence. 

b.  Determine  the  probabilities  of  transition  from  one  state 
to  another  during  periods  of  12  and  24  hours;  i.e. , 
normal  propagation  to  trapping,  super-refraction  to 
normal,  etc. 


I 


Task  Three:  Typical  synoptic  patterns,  frequency,  persistence  and 
transition  probabilities.  Development  of  a capability  to 
assess  refractive  structures  from  synoptic  patterns  in  terms 
of  the  probability  of  the  existing  of  trapping,  super-refractive 
or  sub-refractive  layers,  probable  thickness  of  such  layers, 
the  probable  intensity  of  the  layers  and  a probability  estimate 
(forecast)  of  the  persistence  of  such  layers. 

This  Technical  Report  presents  the  monthly  climatology  of  refractive 
conditions  developed  in  accordance  with  Task  One,  and  measures  of  the 
persistence  and  transition  probabilities  developed  in  accordance  with  Task 
Two. 
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2. 


THE  DATA  BASE 


The  data  available  for  the  development  of  the  required  climatology 
was  in  one  of  two  different  forms: 

a.  CDC  3100  Refractive  Index  Format 

This  stemmed  from  data  held  on  magnetic  tape  by  the  National 
Record  Center  in  the  505  and  645  format,  containing  standard  level  and 
significant  level  information  for  about  15%  of  the  picket  ship  radiosonde 
reports.  This  data  was  processed  by  NEPRF  on  a CDC  3100  computer  to 
obtain  the  refractive  conditions  of  each  layer  and  the  results  recorded  on 
magnetic  tape.  The  format  of  this  refractive  index  information  is  shown 
in  the  Appendix. 

b.  CDC  6500  Refractive  Index  Format 

As  indicated  above,  only  about  15%  of  the  total  picket  ship  data 
was  readily  available  on  magnetic  tape.  The  other  85%  was  in  the  form 
of  plotted  soundings.  These  were  digitized  and  then  processed  or^  a 
CDC  6500  computer  to  obtain  the  refractive  conditions  of  each  layer. 

The  results  were  recorded  on  magnetic  tape  in  the  format  shown  in  the 
Appendix. 

As  part  of  the  overall  data  processing  procedure  the  data  base  was 
converted  to  a common  format.  This  was  combined  with  the  correction  of 
errors  in  the  data  base  and  is  discussed  in  Section  4. 


3. 


THE  REFRACTIVE  PARAMETERS  AND  PROFILE  CLASSIFICATIONS 


3 . 1 Definition  of  Parameters 

The  parameter  N is  simply  a transformation  of  the  variable  speed, 
v , of  the  electromagnetic  propagation  in  the  atmosphere: 


N = 107  ( ~ - 1 


where  the  constant  c is  the  speed  in  a vacuum.  The  formula  applied  to 
calculate  N for  each  reported  level  of  a sounding  takes  the  following 
form: 


N = 776  ~ + 3.73  x 106  ^ 


where  P is  the  pressure  in  mbs,  e is  the  partial  pressure  in  mbs  due  to 
water  vapor,  and  T is  the  temperature  in  degrees  Kelvin. 

N may  be  transformed  into  a parameter  M : 

M = N + 1.57  H (■ 

where  H is  height  above  sea  level  in  meters.  This  new  parameter  merely 
shifts  the  value  of  the  vertical  derivative: 


. dN  + 

dH  dH 


The  variable  speed  of  propagation  can  be  recovered  from  parameter  N: 


1 *-  10"7  N 


The  instantaneous  direction  of  propagation  of  an  infinitesimal 
element  of  an  electromagnetic  wavefront  is  defined  by  the  positive  normal 
to  the  wavefront.  A ray  may  be  defined  as  the  trajectory  of  such  an 
element.  The  trajectory  is  largely  determined  by  the  process  of  refraction. 
The  governing  equation  for  the  tracing  of  a ray  in  a two-dimensional  vertical 
plane  of  the  atmosphere  may  be  expressed  by 

d'7  _ cos  y + I 5v 

ds  ~ R *■  H v >n  W 

where  the  angle  y defines  the  direction  of  propagation  measured  positive 
counter-clockwise  from  the  local  horizontal,  R is  the  earth's  radius  at 
local  sea  level,  s is  the  linear  measure  along  the  trajectory,  and  n is 
the  linear  measure  along  the  wavefront,  with  positive  being  to  the  right 
of  the  trajectory.  The  first  of  the  two  terms  on  the  right-hand  side 
determines  the  apparent  bending  of  the  ray  due  to  the  curvature  of  the 
reference  horizontal.  The  second  term  determines  the  rate  of  refraction. 

The  curvature  term  can  be  eliminated  by  introducing  a modified 
speed  parameter: 


v 


r 


R 

v 

R + H 


(7) 


Substitution  for  v in  Eq.  (6)  and  using  cos  y = - dH/dn  (see  Fig.  1) 
yields: 


dy 

ds 


(8) 
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The  speed  v generally  has  a much  greater  variability  with  height 
in  the  atmosphere  than  it  has  in  any  horizontal  plane.  In  applications 
for  which  the  horizontal  variations  are  negligible,  may  be  treated  as  a 
function  of  height  only.  The  modified  speed  v^  is  then  also  a function 
of  height  only,  and  Eq.  (8)  may  be  transformed  into  a perfect  differential, 
producing  Snell's  law: 


= constant  along  any  one  ray 


This  transformation  into  Snell's  law  is  readily  produced  by  use  of  the 


relationships  shown  in  Fig.  1. 


L\ 


cos  y 


ds  = 


sin  y 


Figure  1.  Relationships  between  the  Local  Differentials 


Thus,  Eq.  (8)  yields 


sin  y 


cos  y 


c 


v 


which  is  the  perfect  differential  along  the  ray 

jL  / cos_i\  = 

“i  \ ) 

All  classifications  of  the  refractive  conditions  of  atmospheric 
profiles  are  based  on  a layered  atmospheric  structure.  Each  reported 
pressure  level  defines  an  interface  between  layers,  and  v^  is 
approximated  to  be  linear  in  H between  evaluations  at  successive 
interfaces . 

That  section  of  a ray  which  falls  into  an  atmospheric  layer  for 
which  dvr/dH  = 0,  appears  to  have  no  bending  in  the  horizontal  frame  of 
reference.  If  dvf/dH  is  positive  or  negative  then  the  ray  appears  to  bend 
downward  or  upward,  respectively,  defining  a circular  arc  in  its  transit  of 
the  layer. 

Differentiation  of  Eq.  (7)  with  height  yields 


dv 

r _ R dv  _ R ^ dv  c 

dH  R + H dH  (r+h)2"  V * dH  ~ R 


Differentiation  of  Eq.  (1)  with  respect  to  height  yields 


dN  . _7  c dv  IQ7  dv 

dH  10  TT.  ~ * T dH 

v dh 


(10) 


(ID 


The  relationship  between  dM/dH  and  dN/dH  is  given  by  Eq.  (4).  Equations 
(4),  (10)  and  (11)  interrelate  the  four  differently-scaled  measures  of  the 
vertical  derivative  of  the  electromagnetic  speed  of  propagation. 


t 
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3.2 


Classification  of  Refractive  Conditions 


3.2.1  Profile  Conditions  Classified  According  to  Layer  Categories 

Profiles  may  be  classified  according  to  the  refractive  categories  of 
the  component  layers . 

The  United  States  Navy  has  separated  the  possible  refractive 
conditions  of  layers  of  the  atmosphere  into  four  discrete  categories  which 
are  termed  trapping  (ducting),  super-refractive,  normal,  and  sub-refractive 
conditions.  For  example,  trapping  conditions  exist  in  a layer  of  the 
atmosphere  if  the  gradient  of  N with  height  is  less  than  -1.57  meters  * or 
equivalently  if 


c dH 

is  greater  than  zero. 

Table  1 shows  the  demarcation  values  for  each  of  the  four  discrete 
categories  in  terms  of  the  four  vertical  derivatives  of  speed. 

The  demarcation  between  Trapping  and  Super-Refraction  corresponds 
to  zero  apparent  refraction  in  the  horizontal  frame  of  reference.  The  other 
two  demarcations  separate  ranges  of  upward  refraction. 


C 


i 
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3.2.2  Profile  Conditions  Classified  by  Continuous  Parameters 

A technique  for  representing  the  refractive  conditions  of  a profile 
is  to  separate  it  into  layers,  the  vertical  extent  of  each  layer  being 
determined  by  the  appropriate  demarcation  values  given  in  Table  1. 

Thus  each  layer  is  assigned  to  one  of  the  four  discrete  categories. 
Although  such  a discontinuous  representation  can  be  detailed  and  is 
operationally  significant,  it  nevertheless  lacks  the  continuity  desirable 
for  some  aspects  of  our  investigations. 

We  therefore  have  designed  two  continuous  parameters  as 
supplementary  definitions  of  the  refractive  conditions  of  a profile. 
Parameter  S is  a measure  of  the  strength  of  surface  ducting  or  of  the 
proximity  to  surface  ducting.  Parameter  D is  a measure  of  the  strength 
of  total  tropospheric  ducting  or  of  the  proximity  to  ducting. 

The  definition  of  the  parameter  S is 


S = 


7 

— r~  ( v . - v ) provided  v . > v 

; H*  \ r, i r,o ) r,i  r,o 


where  vr  Q is  the  surface  value  of  the  modified  speed  profile,  vr  , and 
vf  i is  the  maximum  value  in  the  profile  above  the  surface.  H*  is  a 
specified  scaling  height.  The  extension  of  the  definition  in  the  case  of 
no  surface  duct,  i.e. , if  vr  Q is  the  maximum  speed  in  the  profile,  is 


S * 


where 


and  Hi  is  the  height  of  level  i . This  extension  gives  a measure  of  the 
proximity  of  the  profile  to  a surface  duct  situation.  This  negative  extension 


of  S is  not  necessarily  related  to  the  bottom  layer  defined  by  the  first 
reported  level  (i  = 1)  above  the  surface;  the  maximum  may  occur  for  i = 2 
or  higher.  The  interpretation  of  a negative  value  of  S is  the  same  as  for 
the  first  row  of  demarcations  given  in  Table  1. 

A positive  value  of  S is  interpreted  by  application  of  Snell’s  law 
for  the  limiting  ray  at  the  surface: 


(Limit)  cos  V 


V . - V V . - V 

= 1 _ LlO  w _ _Ld 

V . c 

r,  i 


*1-10  H^S 


The  larger  is  S the  greater  is  the  limiting  angle  for  surface  ducting. 

The  total  tropospheric  ducting  parameter  D is  similar  to  S except 
that  it  need  not  be  surface  based: 


_ i°7  i / \ 

D = ~ I v , ~ v , 

c H*  \ rj  r.i ) 


provided  v > v . and  j > i 
r,J  r,i 


D » 

c \ dH  / 


for  no  ducting. 


■ 

l 


The  interpretation  of  D for  no  ducting,  that  is  for  D negative, 
is  the  same  as  the  first  row  of  Table  1.  A positive  value  of  D is 
interpreted  by  application  of  Snell's  law  for  the  limiting  ray  at  level  i 


(Limit)  cos  y 


For  the  purpose  of  presenting  the  statistics  of  persistence*  and 
transition,  it  is  convenient  to  establish  discrete  categories  for  the  S and 
D parameters.  These  are  shown  in  Table  2,  together  with  the  first  row  of 
Table  1 for  comparison  purposes  . 


Persistence  probabilities  are  a sub-set  of  transition  probabilities;  i.e. , 
that  sub-set  where  no  transition  from  one  refractive  category  to  another 
occurs  in  the  space-time  envelope  under  consideration. 


Table  2 


Discrete  Categories  of  the  S and  D Parameters 


Trapping/Ducting 

Super-Refraction 

Normal 

Strong  (DS)  Moderate  (DM)  Weak  (DW) 

Strong  (SS)  Weak  (SW) 

(N) 

,7  dv 
10  r 

c dH 


-0 . 79-**-  -1.57 


-0.40-*-*. -0.79. 


-0. 40-*-* -0.79. 


The  S and  D parameters,  calculated  in  accordance  with  Eqs.  (12), 
(13),  (15)  and  (16),  are  a measure  of  the  most  positive/least  negative  values. 
As  a consequence,  no  values  of  S and  D corresponding  to  sub- refractive 
conditions  were  ever  observed. 

Every  profile  provides  a value  for  S and  a value  for  D . Because 
D is  determined  from  the  whole  profile,  whereas  S is  only  applicable  to 
the  surface,  the  values  of  S and  D will  coincide  if  the  base  of  the  most 
positive/least  negative  layer  is  at  the  earth's  surface.  However,  if  this 
layer  is  aloft,  then  the  profile  will  give  different  values  for  S and  D. 


For  ducting/trapping  conditions  (eqs.  (12),  (15) ) the  parameters 
S and  D are  directly  related  to  the  limiting  rays  (Eqs.  (14),  (17)).  The 
quantity  H*  is  included  in  the  definition  of  positive  S and  D values,  not 
in  order  to  express  them  as  gradients,  but  to  express  them  in  a scale 
similar  to  their  negative  definitions. 
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4. 


PROCESSING  THE  DATA  BASE 


After  construction  of  a data  listing,  the  position  given  with  each 
report  was  plotted;  those  reports  which  either  fell  over  land,  were  too 
close  to  land,  or  which  were  geographically  remote  from  other  reports, 
were  considered  suspect.  Such  reports  were  checked  by  a variety  of 
methods  including  reference  to  the  original  plotted  soundings.  A search  was 
made  for  duplicate  reports  for  subsequent  elimination.  The  maximum  and 
minimum  base  levels  and  thicknesses  of  refractive  layers  were  extracted 
and  suspect  extreme  values  checked.  Finally  a cross-check  was  made 
of  the  positions  and  date-times  of  the  reports  on  the  plotted  soundings 
against  the  same  information  held  on  magnetic  tape  after  digitization. 

The  errors  revealed  by  this  search,  and  the  associated  corrective 
procedures,  are  tabulated  in  Table  3. 
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4.2  Data  Base  Reformatting  and  Sorting 

As  explained  in  Section  2,  the  basic  data  was  held  in  one  of  two 
different  formats.  A third  format,  more  appropriate  to  the  investigation, 
was  devised  and  conversion  programs  were  written  and  run  for  each  data 
set.  The  final  format  of  the  data  base  is  shown  in  the  Appendix.  The 
reformatting  programs  corrected  the  errors  in  the  data  base  revealed  by 
the  error  search  (see  Section  4.1)  and  also  computed  and  incorporated 
the  S and  D parameters  (see  Section  3.2.2). 

The  reformatted  data  was  then  sorted  and  merged  into  a single  data 
base,  including  the  S and  D parameters.  The  basic  ordering  of  the  data 
was  by  time,  by  latitude  and  then  by  longitude.  The  merged  and  sorted 
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Table  3 


Data  Base  Errors  and  Corrections 


Corrections  to  3100  Tape 


Original  Number 
Corrections  to  Latitude 
Corrections  to  Data 

Bad  Data — Eliminated 
Duplicates — Eliminated 
Number  Eliminated 
Number  Retained 


985 


1 

15 


976 


Corrections  to  6500  Tape 


Original  Number 

Corrections  to  Latitude/Longitude  68 

Corrections  to  Date  138 

Corrections  to  Ship  Number  9 

Corrections  to  Data  10 

Bad  Reports  (No  matching 
Plotted  Soundings) — Eliminated 
Bad  First  Levels — Eliminated 
Duplicates  — Eliminated 
Number  Eliminated 
Number  Retained 

Total  Number  Retained 


5987 


10 

16 

15 


41 


5946 

6922 
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data  base  contained  6922  profiles  (00Z  or  12Z  only)  on  magnetic  tape  in 
the  format  shown  in  the  Appendix.  These  profiles  provided  the  input  to  the 
statistical  analysis  described  in  Section  5. 


4 . 3 The  Time:  Latitude  Distribution  of  Profiles 

Figures  2 through  6 show  the  chronological  distribution  of  profiles 
in  the  data  base  as  a function  of  latitude;  each  figure  shows  the  data  for 
two  years.  Time  is  shown  along  the  abscissa,  the  numbers  above  the 
months  being  the  number  of  profiles  for  that  month.  Latitude  is  shown 
along  the  ordinate  to  the  nearest  0.1°  of  latitude.  Individual  reports  are 
indicated  by  a short  vertical  line  centered  on  the  time-latitude  coordinates 
of  the  report.  The  numbers  between  whole  latitude  lines  are  the  number 
of  profiles  for  that  degree  of  latitude  for  the  two-year  period.  The 
numbers  appearing  between  the  latitude  lines  on  the  right-hand  side  of 
Fig.  6 give  the  total  number  of  profiles  occurring  in  that  latitude  band  for 
the  10-year  period  of  the  data  base. 
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5. 


STATISTICAL  COMPILATIONS 


5 . 1 Establishment  of  Latitude  Zones  for  the  Statistical  Analysis 

In  order  to  detect  latitudinal  variations  in  refraction  conditions,  it 
was  necessary  to  separate  the  profiles  in  the  data  base  into  latitude  zones. 
The  demarcation  latitudes  selected  (29.0°N,  34.5°N,  37.5°N,  41.0°N, 
46.03N,  51.0°N)  were  a compromise  involving  several  considerations  such 
as: 

a.  having  a reasonable  number  of  profiles  in  each  zone,  these 
profiles  also  being  representative  in  time 

b.  having  roughly  equal  zone  sizes; 

c.  minimizing  the  number  of  reports  actually  on  the  selected 
lines  of  demarcation; 

d.  minimizing  the  movement  of  ships  (making  regular  reports) 
across  the  selected  lines  of  demarcation; 

e.  avoiding  concentrations  of  data  near  the  zone  boundaries. 

The  compromise  was  effected  by  plotting  the  number  of  reports 
occurring  in  0.5°  latitude  x 0.5  longitude  areas  for  each  month  and  for 
all  months  combined,  and  considering  these  plots  (Figs.  7 through  19)  in 
conjunction  with  Figs.  2 through  6.  On  Figs.  7 through  19,  the  finally- 
selected  demarcation  latitudes  are  shown  as  dashed  lines;  the  numbers 
along  the  ordinate  between  the  latitude  digits  are  the  number  of  reports 
in  that  zone  for  that  month,  and  the  number  at  the  bottom- left  of  the  figures 
is  the  total  number  of  reports  for  all  zones. 
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5.2 


Statistical  Development  of  the  Monthly  Climatology  of  Refractive 
Conditions  (Task  One) 


Because  of  the  operational  significance  of  surface  ducting,  surface 
and  elevated  ducts  have  been  considered  as  separate  categories.  The 
following  table  shows  the  parameters  used  to  specify  each  significant 
category  of  refraction: 


Table  4 

Specifying  Parameters  for  the  Discrete  Categories 
of  Refraction 


Non-Surface  Ducts 

Trap  Base  Height 
Trap  Thickness 
Trap  Gradient 
Duct  Thickness 
Duct  Base 


Surface  Ducts 

Trap  Base  Height 
Trap  Thickness 
Trap  Gradient 
Duct  Thickness 


Super-Refractive 
Layers 

Base  Height 

Thickness 

Gradient 


Sub-Refractive 

Layers 

Base  Height 
Thickness 


i 

f 


l 
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The  following  definitions  are  relevant  to  the  specifying  parameters: 


! 


Trapping  Layer: 


A vertical  extent  in  which  dM/dH  is  negative 
throughout.  The  full  extent  of  such  negative 
gradient  defines  one  trapping  layer.  A single 
profile  may  include  several  trapping  layers , 
separated  by  regions  in  which  dM/dH  is  positive 
throughout. 


Trap  Base  Height:  The  height  of  the  bottom  of  the  trapping  layer. 


L 


i 
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Height  of  the  top  minus  height  of  the  bottom  of 
a trapping  layer. 


Trap  Thickness: 


Trap  Gradient: 
Duct  Base: 


Duct  Thickness: 


The  mean  gradient,  dM/dH , of  a trapping  layer. 

One  duct  base  corresponds  to  each  trapping  layer. 
For  a non-surface-based  duct  the  base  is  at  the 
first  level  encountered,  in  progressing  downward 
from  the  base  of  the  trapping  layer,  where  the 
speed  equals  the  speed  at  the  top  of  the  trapping 
layer.  If  the  surface  is  encountered  before  this 
condition  is  met  then  the  duct  is  surface  based. 

Height  of  the  top  of  a trapping  layer  minus  the 
height  of  its  duct  base.  In  the  case  of  a surface- 
based  duct,  this  thickness  is  the  height  of  the  top 
of  the  associated  trapping  layer. 


Super-Refractive  Layer:  A vertical  extent  in  which  dM/dH  lies  between 

0 and  0.79  throughout.  A single  profile  may 
include  several  super-refractive  layers,  each 
having  a base  height  (the  lower  extent  of  the  layer) 
a thickness,  and  a mean  gradient. 

Sub-Refractive  Layer:  A vertical  extent  in  which  dM/dH  is  greater  than 

1.57  throughout.  A single  profile  may  include 
several  sub-refractive  layers,  each  having  a base 
height  and  a thickness. 


Table  5 summarizes  the  parameters  and  statistical  measures  used  to 
produce  the  climatologies  of  refractive  conditions  developed  in  accordance 
with  the  requirements  of  Task  One. 


1 


C 


-20- 


VN, 


< II 


*■ 


k 


Table  5 

Parameters  and  Statistical  Measures  Used  in 
Developing  the  Monthly  Climatologies  of  Refractive  Conditions 


Statistical  Measures  mean,  standard  deviation, 

percentage  frequency  of  occurrence  3 

Specifying  Parameters  see  Table  4 14 

Months  12  calendar  months  plus 

1 all-months  category  13 

Times  00Z  and  12Z,  plus  00Z  & 12Z 

combined  3 

Latitude  Zones  5 zones  plus  1 all-zones  category  6 


The  various  charts  and  tables  used  to  present  the  monthly  climatologies  of 
refractive  conditions  are  described  in  Section  5.4. 

Note: 

a.  With  reierence  to  Table  4,  a single  profile  may  contain  any 
number  (within  reason)  of  surface  ducts,  non-surface  ducts,  super-refractive 
layers  and  sub-refractive  layers,  and  these  may  all  co-exist  in  a given 
profile.* 

b.  With  reference  to  Table  5,  the  statistics  for  00Z  and  12Z  were 
produced  separately  and  in  combination  to  detect  diurnal  biases. 

c.  The  basic  flow  chart  of  the  computer  program  used  to  produce 
the  monthly  climatologies  of  refractive  conditions  is  given  in  the  Appendix. 


I 
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See  Table  20. 
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5 . 3 Statistical  Development  of  Persistence  and  Transition  Probabilities 

for  Refractive  Conditions  (Task  Two) 

The  statistics  for  the  probability  of  transition,  from  one  category 
of  refraction  to  another,  including  the  percentage  frequency  of  occurrence, 
were  developed  using  the  S and  D parameters  defined  and  discussed  in 
Section  3.2.2.  Transition  (or  persistence)  probabilities  between  the  six 
discrete  categories  of  refraction  for  these  parameters  laid  down  in  Table  2 
were  determined  by  considering  changes  either  in  space  or  in  time,  but  not 
space-and-time.  For  detecting  any  seasonal  variations  the  year  was 
divided  into  four  quarters  designated  Winter,  Spring,  Summer  and  Fall,  but 
no  true  seasonal  regime  changes  should  be  inferred. 

Spatial  transition  probabilities  were  determined  between  adjacent 
latitude  zones  for  a constant  time;  temporal  transition  probabilities  were 
determined  in  each  latitude  zone,  including  one  all-zones  category,  for 
12  hours  and  24  hours  in  the  future  based  on  starting  times  of  00Z  and  12Z. 

Table  6 summarizes  the  parameters  and  statistical  measures  used 
in  developing  the  required  transition/persistence  probabilities.  Specific 
details  of  the  method  used  for  the  compilation  of  each  set  of  statistics  are 
given  in  Section  5.4. 


Table  6 

Parameters  and  Statistical  Measures  Used  in  Developing 
the  Statistics  of  Transition  and  Persistence 


Statistical  Measures  number  of  profiles,  percentage  frequency 

of  occurrence,  cumulative  percentage 
frequency  of  occurrence 

Specifying  Parameters  S and  D , each  with  6 discrete  categories 

of  refraction;  i.e.,  DS,  DM,  DW,  SS, 

SW  and  N.  (See  Section  3.2.2  and 
Table  2) 


Period 


4 quarters  of  the  year: 

January  through  March  ("Winter") 
April  through  June  ("Spring") 

July  through  September  ("Summer") 
October  through  December  ("Fall") 


Times  (for  temporal  transition  probabilities) 


00Z  compared  with  00Z  +12  hours 
and  OOZ  +24  hours 
12Z  compared  with  12Z  +12  hours 
and  12Z  +24  hours 


(Temporal  transition  probabilities  were  determined  within  each  of 
the  5 latitude  zones,  plus  1 all-zones  category.) 

Latitude  Zones  (for  spatial  transition  probabilities) 

5 latitude  zones , each  compared  with 
the  latitude  zone  to  the  north  (4)  and 
to  the  south  (4) 

(Spatial  transition  probabilities  were  determined  for  simultaneous 
occurrences  only.) 
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Description  of  Tables  and  Figures 


This  section  described  the  format  of  the  tables  and  figures  used  to 
present  the  statistical  summaries  of  refractive  conditions. 

5.4.1  Monthly  Mean  Values  for  Non-Surface  Ducts 

Figures  20  through  24  show  the  mean  values  of  the  five  specifying 
parameters  for  non-surface  ducts  (see  Table  4).  Each  figure  contains  six 
histograms,  one  for  each  latitudinal  zone  and  one  for  all-zones  combined  — 
the  top-right  histogram  in  each  figure.  The  relevant  latitudes  are  written 
along  the  right-hand  edge  of  each  histogram.  The  ordinate  shows  either 
height,  thickness  or  gradient;  the  units  of  height  and  thickness  are  meters, 

and  the  units  of  gradient  are  M units  per  100  meters.  Mean  values  are  * 

( 

given  for  00Z,  12Z,  and  combined  00Z  & 12Z.  NS  refers  to  the  number  of 
available  soundings  and  is  given  for  each  calendar  month  and  all-months. 

N refers  to  the  number  of  occurrences  of  the  specifying  parameter  for  each 
time  and  latitude  zone.  Note  that  more  than  one  of  these  occurrences  may 
have  been  contributed  by  an  individual  sounding  (see  Section  5.2). 

5.4.2  Monthly  Mean  Values  for  Surface  Ducts 

Figures  25  through  28  show  the  mean  values  of  the  four  specifying 
parameters  for  surface  ducts  (see  Table  4).  See  Section  5.4.1  for 
description  of  contents  of  these  four  figures. 

5.4.3  Monthly  Mean  Values  for  Super-Refractive  Conditions 


Figures  29  through  31  show  the  mean  values  of  the  three  specifying 
parameters  for  super-refractive  layers  (see  Table  4).  See  Section  5.4.1 
for  description  of  contents  of  these  three  figures. 
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5.4.4  Monthly  Mean  Values  for  Sub-Refractive  Conditions 

Figures  32  and  33  show  the  mean  values  of  the  two  specifying 
parameters  for  sub-refractive  layers  (see  Table  4).  See  Section  5.4.1  for 
description  of  contents  of  these  two  figures. 


5.4.5. 1 Percentage  Frequency  Distribution  of  Trap  Base  Heights 

Figures  34  through  39  show  the  percentage  frequency  distribution  of 
trap  base  heights  for  all  occurrences  of  trapping,  both  non-surface  and 
surface  based.  Each  figure,  one  for  each  latitude  zone,  and  one  for  all-zones 
combined,  shows  thirteen  histograms,  one  for  each  calendar  month  and  one 
for  all-months  combined.  In  each  histogram  the  ordinate  is  the  base  height 
of  the  trap  in  meters;  the  abscissa  is  the  percentage  frequency.  The  latter 
was  evaluated  for  each  height  interval  as  a percentage  of  t.he  total 
occurrences  of  trapping  for  that  particular  month  and  latitude  zone.  Thus, 
for  example,  referring  to  the  June  histogram  in  Fig.  34,  of  all  the  recorded 
occurrences  of  trapping  layers  in  June  1957-1965  between  the  latitudes  of 
29.0°  and  34.5°N,  10%  of  them  had  base  heights  of  250-500  meters. 

5. 4. 5. 2 Percentage  Frequency  Distribution  of  Trap  Thicknesses 

Figures  40  through  45  are  basically  similar  to  those  described  in 
Section  5.4.5. 1,  the  difference  being  that  each  ordinate  shows  the  percentage 
frequency  distribution  of  the  thicknesses  of  the  traps  rather  than  the  base 
heights . 

5.4.6. 1 Percentage  Frequency  Distribution  of  Duct  Base  Heights 

Figures  46  through  51  are  similar  to  those  described  in  Section  5.4.5. 1 
but  show  the  percentage  frequency  distribution  of  the  heights  of  the  duct 
base. 
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5. 4. 6. 2  Percentage  Frequency  Distribution  of  Duct  Thicknesses 

Figures  52  through  57  are  similar  to  those  described  in  Section 

5.4.5. 1  but  each  ordinate  shows  the  percentage  frequency  distribution  of 
the  thicknesses  of  the  ducts. 

5.4.7. 1 Percentage  Frequency  Distribution  of  the  Base  Heights 
of  Super-Refractive  Layers 

Figures  58  through  63  are  similar  to  those  described  in  Section 

5.4.5. 1  but  show  the  percentage  frequency  distribution  of  the  base  heights 
of  super-refractive  layers. 

5.4. 7. 2 Percentage  Frequency  Distribution  of  the  Thicknesses 
of  Super-Refractive  Lavers 

Figures  64  through  69  are  similar  to  those  described  in  Section 

5.4.5. 1  but  show  the  percentage  frequency  distribution  of  the  thicknesses 
of  super-refractive  layers. 

5.4.8. 1 Percentage  Frequency  Distribution  of  the  Base  Heights 
of  Sub-Refractive  Lavers 

Figures  70  through  75  are  similar  to  those  described  in  Section 

5.4.5. 1  but  show  the  percentage  frequency  distribution  of  the  base  heights 
of  sub-refractive  layers. 

5 . 4 . 8 . 2 Percentage  Frequency  Distribution  of  the  Thicknesses 
of  Sub-Refractive  Lavers 

Figures  76  through  81  are  similar  to  those  described  in  Section 

5.4.5. 1  but  show  the  percentage  frequency  distribution  of  the  thicknesses  of 
sub- refractive  layers. 
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5.4.9  Statistical  Summaries  for  the  Categories  of  Refraction 

5.4.9. 1 Tables  7 through  19  show  the  means  and  standard  deviations 
of  the  fourteen  specifying  parameters  detailed  in  Table  4. 

As  an  example,  consider  Table  7.  This  table  is  for  all  January 
observations  and,  for  each  latitude  zone  and  for  all-zones  combined,  the 
specifying  parameters  for  each  discrete  category  of  refraction  are  listed  in 
the  left-hand  column.  The  first  three  quantities  (under  "NUMBER")  show 
the  number  of  occurrences  of  the  four  discrete  categories  of  refraction.  For 
example,  considering  the  zone  29.0°N  to  34. 5° N,  there  were  14  super- 
refractive  layers  occurring  in  all  Januaries  at  00Z  and  19  at  12Z,  giving  a 
total  of  33.  Any  particular  profile  may  have  had  more  than  one  super- 
refractive  layer,  and  the  next  three  quantities  (under  "NUMB  RPTS  1 OR 
MORE")  show  the  number  of  profiles  exhibiting  super-refractive  layer(s). 

Thus,  in  this  example,  there  were  11  profiles  with  one  or  more  super- 
refractive  layer(s)  in  all  Januaries  at  00Z  and  14  at  12Z,  giving  a total 
of  25. 

For  this  particular  sub-table  the  top  line  shows  there  were  24  reports 
at  00Z  and  25  at  12Z,  a total  of  49.  Under  "PERCENT  RPTS  1 OR  MORE"  the 
number  of  profiles  exhibiting  super-refractive  layers  are  expressed  as 
percentages  of  the  appropriate  number  of  profiles  available.  Thus  45%  of 
January  profiles  in  this  latitude  zone  had  one  or  more  super-refractive 
layers  at  00Z,  56%  at  12Z,  and  51%  at  either  00Z  or  12Z. 

The  three  quantities  under  "MEANS"  show  the  r san  values  of  the 
specifying  parameters  at  00Z,  12Z,  and  00Z  & 12Z  combined,  the 
appropriate  standard  deviations  being  shown  in  the  final  three  columns. 

The  means  were  calculated  for  all  occurrences.  Thus,  using  the  same 
sub-table,  the  mean  height  of  all  33  super-refractive  layers  occurring  in 
January  profiles  for  the  latitude  zone  29.0°N  to  34. 5° N was  1162  meters 
with  a standard  deviation  of  798  meters. 
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Note:  In  Tables  7 through  19  , the  term  "Intensity"  is  synonymous  with 
"Mean  Gradient";  for  trapping  conditions  the  sign  of  this  gradient  has  been 
reversed.  The  units  of  intensity  are  M units  per  100  meters. 

5. 4. 9. 2 Table  20  shows  the  numbers  of  occurrences  of  profiles 
exhibiting  up  to  six  (for  each  category)  co-existing  ducts,  super-refractive 
and/or  sub-refractive  layers. 

5.4.10  Seasonal  Distribution  of  Parameter  S and  Parameter  D 

Tables  21  through  24  show  the  seasonal  distributions  (see  Table  6) 
of  the  S and  D parameters  for  each  of  the  six  latitude  zones,  00Z  and  12Z 
reports  being  shown  separately.  Values  of  S and  D are  shown  in  terms  of 
the  discrete  ranges  laid  down  in  Table  2. 

Thus,  for  example.  Table  21  shows  the  distribution  of  Parameter  S 
derived  from  00Z  reports  (profiles).  Considering  the  two  sub-tables  relating 
to  the  latitude  zone  29.0°N-34.5°N,  the  left-hand  sub-table  shows  that 
in  Spring  (April  through  June) , there  were  170  relevant  profiles  in  the  data 
base.  Separation  of  these  170  profiles  into  the  discrete  ranges  of  refraction 
gave  93  cases  of  normal  (N)  refraction,  39  cases  of  weak  super-refraction 
(SW) , and  so  on.  The  right-hand  sub-table  shows  these  profiles  expressed 
as  a (truncated)  percentage  frequency  of  occurrence. 

5.4.11  Comparison  of  the  Distributions  of  Parameters  S and  D in 
Adjacent  Latitude  Zones 

Tables  25  through  28  compare  the  distributions  of  Parameters  S and 
D in  one  latitude  zone  with  the  corresponding  distribution  in  an  adjacent 
latitude  zone. 

The  sub-set  of  the  data  base  required  for  the  comparison  was 
determined  by  considering  two  adjacent  latitude  zones  and  time-matching 
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the  profiles  from  each  zone.  Pairs  of  simultaneously-occurring  profiles 
obtained  in  this  way  were  used  to  form  the  distributions  of  Parameters  S 
and  D in  the  two  latitude  zones . Insufficient  matched  pairs  were  obtained 
in  each  season  to  provide  a statistically-meaningful  sample,  and  therefore 
the  comparative  distribution  of  S and  D were  derived  for  all-seasons 
combined. 

As  an  example,  the  top-left  sub-table  in  Table  25  shows  latitude 
zone  29.0°N-34.5°N  compared  with  latitude  zone  34.5°N-37.5°N  for  the 
Parameter  S.  Addition  of  the  right-hand  column  of  the  sub-table  (under 
"TOT")  shows  that  there  were  451  corresponding  pairs  of  simultaneously 
observed  profiles.  Of  these  451  pairs  of  soundings,  230  profiles  fell  in 
the  "N"  category  (see  Table  2)  in  latitude  zone  29.0°N-34.5°N.  The  230 
corresponding  profiles  in  latitude  zone  34 . 5°  N-37. 5° N had  the  following 
distribution:  N — 175  profiles;  SW — 32  profiles;  SS — 3 profiles;  DW — 17 
profiles;  DM — 3 profiles;  and  SS — 0 profiles.  The  right-hand  sub-table 
shows  these  profiles  expressed  as  a (truncated)  percentage  frequency  of 
occurrence . 

Considering  now  Table  26,  the  top  left-hand  sub-table  shows  the 
same  451  corresponding  pairs  of  simultaneously  observed  profiles.  Of 
these  451  pairs  of  soundings,  322  fell  into  category  N in  latitude  zone 
34. 5° N-37. 5° N.  The  322  corresponding  profiles  in  latitude  zone  29.0°N- 
34.5°N  had  the  following  distribution:  N — 175  profiles;  SW — 83  profiles; 
SS — 23  profiles;  DW — 21  profiles;  DM — 10  profiles;  and  DS — 10  profiles. 
The  right-hand  sub-table  shows  these  profiles  expressed  as  a percentage 
frequency  of  occurrence. 

The  remainder  of  the  sub-tables  in  Tables  25  and  26  are  similar, 
comparing  the  distributions  of  Parameter  S in  other  adjacent  latitude 
zones.  Tables  2 7 and  28  are  similar  to  Tables  25  and  26,  but  show  the 
distribution  of  Parameter  D . 
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Figure  82  is  a scattergram  for  the  same  451  pairs  of  profiles 
discussed  in  the  above  example.  The  ordinate  and  abscissa  show 
Parameter  S as  a continuous  variable  for  the  two  adjacent  latitude  zones, 
the  S-values  of  corresponding  pairs  of  profiles  being  plotted  on  the 
scattergram.  Similar  figures  were  plotted  for  Parameters  S and  D for  all 
adjacent  latitude  zones;  they  are  not  included  in  this  Report  as  the 
tabular  format  is  more  informative . 


5.4.12  Transition  and  Persistence  Probabilities  for  Parameters 
S and  D with  Time 

Tables  29  through  61  show  how  the  observed  distributions  of 
Parameters  S and  D changed  (or  persisted)  in  time  within  the  same  latitude 
zone.  Tables  29  through  44  refer  to  Parameter  S , while  Tables  45  through 
60  refer  to  Parameter  D . 

The  sixteen  tables  for  Parameters  S are  divided  into  four  sets — one 
set  for  each  of  the  four  "seasons"  (see  Table  6).  The  four  tables  in  each  of 
these  sets  show  the  changes  in  distribution  of  Parameter  S after  a time-lapse 
of  12  hours  based  on  00Z  and  on  12Z  soundings,  and  after  a time-lapse  of 
24  hours  based  on  00Z  and  on  12Z  soundings.  The  sixteen  tables  for 
Parameter  D are  organized  similarly. 

The  sub-sets  of  the  data  base  required  to  determine  these  changes 
in  distribution  were  obtained  by  matching  profiles  from  the  base-time  with 
profiles  measured  12  or  24  hours  later  in  the  same  latitude  zone.* 

As  an  example  of  the  contents  of  each  table,  consider  Table  51. 

This  shows  the  change  in  distribution  of  Parameter  D over  a 24-hour  period 
based  on  00Z  soundings  for  April  through  June.  The  change  in  distribution 

Almost  invariably,  these  profile-pairs  were  from  soundings  made  by  the 
same  picket  ship. 


of  Parameter  D is  shown  for  the  five  latitude  zones  and  one  all-zones 
combined  category.  The  base-time  of  the  table  (in  this  case  00Z)  is  shown 
along  the  left-hand  edge  of  the  three  sub-tables  given  for  each  latitude 
zone. 

The  bottom-left  sub-table  in  Table  51  refers  to  the  latitude  zone 
29 . 0° N-5 1 . 0° N;  i.e.,  the  all-zones-combined  category.  Summing  the 
number  of  profiles  under  "TOT"  shows  there  were  788  profiles  in  this 
particular  sub-set  of  the  data  base.  From  the  values  of  Parameter  D 
computed  from  these  788  profiles,  227  (under  "TOT")  fell  into  the  N 
category  (see  Table  2)  at  00Z,  the  base-time  of  the  table.  24-hours  later, 
the  corresponding  profile  distribution  was:  N — 105;  SW — 61;  SS-36; 

DW — 23;  DM — 2;  and  DS — 0.  The  other  rows  of  this  sub-table  show  the 
changes  in  distribution  for  the  other  refractive  categories . 

The  left-hand  sub-tables  show  distributions  based  on  numbers  of 
profiles.  In  the  center  sub-table,  these  numbers  have  been  converted  to 
(truncated)  percentage  frequencies  ol  occurrence.  The  right-hand  sub-table 
of  each  set  of  three  shows  cumulative  percentage  frequencies  of  occurrence, 
summed  from  the  right.  For  example,  the  bottom-right  sub-table  of  Table  51 
shows  that  for  all  profiles  which  were  N at  00Z,  24-hours  later  100%  of 
the  corresponding  profiles  were  N or  stronger,  53%  were  SW  or  stronger, 

2 6%  were  SS  or  stronger,  11%  were  DW,  and  there  were  less  than  1%  in  the 
DM  and  DS  categories.  From  the  same  sub-table  it  can  be  seen  that,  after 
24  hours,  there  is  only  an  11%  chance  of  initially  normal  conditions  turning 
into  a ducting  situation,  whereas  if  there  is  initially  a strong  duct,  there 
is  a 75%  chance  that  ducting  will  still  be  observed  after  24  hours. 

The  contents  of  the  other  tables  are  similar  to  Table  51. 

Figures  83  and  84  show  sample  scattergrams  for  Parameters  S and  D 
for  the  latitude  zone  29 . 0° N-34 . 5°N,  all  seasons  combined,  for  a time- 
lapse  of  24  hours  based  on  00Z  profiles.  Other  similar  figures  were  plotted 
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but  have  not  been  included  in  this  Report  as  the  tabular  format  is  more 


informative 


Various  investigations  and  evaluations  of  practical  use  can  be  based 
on  the  statistical  compilations  and  associated  presentations  contained  in 
this  Report.  The  "pay-off"  of  such  evaluations  will  depend  on  the 
operational  context  and,  to  be  of  maximum  benefit,  should  involve  the 
techniques  of  operational  analysis.  The  prospective  users  of  the  type  of 
information  contained  in  this  Report  therefore  must  themselves  provide  an 
input  to  the  evaluation  process.  As  the  operational  significances  of  the 
variabilities  are  appreciated,  it  will  then  be  possible  to  modify  and 
condense  the  presentation  formats  to  show  more  directly  the  desired 
practical  relationships  between  the  various  parameters  determining  and 
affecting  refractive  conditions. 


Chronological  Distribution  of  Radiosonde  Profiles  with  Latitude 


1958  1959 

Miwro/ya-  aiqq 

Figure  3.  Chronological  Distribution  of  Radiosonde  Profiles  with  Latitude:  1958-1959. 


UHtM 


LATITUDE 


LONGITUDE 


' V i; 


C'  V3 
— *>  \ \ 

2r) 

V,  <r 


1 1 I 

I 7 11 

Wa~ 


I 1|  2 3 


3 23  3 I 3 I 


I 2 11  I 23 
z 13  l 


3 1 13  215  2 


30; 


NUfBFR-ossa 


JANUARY 


Figure  7.  Geographical  Distribution  of  Radiosonde  Reports:  January  1957-1965. 
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Figure  8.  Geographical  Distribution  of  Radiosonde  Reports:  February  1957-1965. 
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Figure  10.  Geographical  Distribution  of  Radiosonde  Reports:  April  1957-1965. 
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Figure  13.  Geographical  Distribution  of  Radiosonde  Reports:  July  1957-1964 
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Figure  14.  Geographical  Distribution  of  Radiosonde  Reports:  August  1957-1964 
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Figure  15.  Geographical  Distribution  of  Radiosonde  Reports:  September  1957-1964. 
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r’igure  18.  Geographical  Distribution  of  Radiosonde  Reports:  December  1956-1964. 
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Figure  19.  Geographical  Distribution  of  Radiosonde  Reports 
All  Months  October  1956-June  1965 


Monthly  Mean  Values  of  Trap  Thickness  for  Surface  Ducts 


Figure  27.  Monthly  Mean  Values  of  Trap  Gradient  for  Surface  Ducts  (sign  of  gradient  reversed) 
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Percentage  Frequency  Distribution  ns=nit3er  cf  scucincs  fdr  tkht  kcnth 

of  Trap  Thicknesses 

Latitude  Zone  4 1 . 0°  N-46 . 0°  N 


Percentage  Frequency  Distribution  NS™*  0F  S0l^IfCS  ™ ™T  "*TH 
of  Duct  Base  Heights 
Latitude  Zone  34.5°N-37.5°N 
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Figure  51.  N =nuhber  of  the  conditions  for  twit  month 

Percentage  Frequency  Distribution  ns=num3er  of  soundings  for  thht  month 

of  Duct  Base  Heights 
Latitude  Zone  29.0  N-51.0  N 


Percentage  Frequency  Distribution  ns=nukser  of  soundings  for  thrt  mchth 
of  Duct  Thicknesses 
Latitude  Zone  37 . 5°  N-4 1 . 0°  N 


Percentage  Frequency  Distribution  ns=nu«3Fr  of  soundings  for  thrt  month 
of  Duct  Thicknesses 
Latitude  Zone  46.0°N-51.0°N 
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Figure  69 . N juicer  of  the  coitions  for  tutt  fcsth 

Percentage  Frequency  Distribution  ns^nuhser  of  soundings  for  thpt  hcntm 

of  the  Thicknesses  of  Super- 

Refractive  Layers 

Latitude  Zone  29 . 0°  N-51 . 0°  N 


450  • 

3ZI  ■ 


lit 


SS3NM3IH1 


SS3N)0JHi 


SS3NM31IU 


110 


N^Z233  2,(6=0165.,  „ j,N=Z371  - (6=0163 „ _ Si?  .,N=B267  a (6=0197 


Figure  81.  h muxazn  of  the  cohutions  for  t«u  month 

Percentage  Frequency  Distribution  ns^nuhser  of  soinjincs  for  thrt  hwth 
of  the  Thicknesses  of  Sub- 
Pefractive  Layers 
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Scattergram  Comparing  the  Distributions  of  Parameter  S in 
Two  Adjacent  Latitude  Zones  for  Simultaneously-Occurring 
Pairs  of  Profiles 
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Figure  83.  Scattergram  Showing  the  Distribution  of  Parameter  S after 
24  Hours  Based  on  00Z  Observations 
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Statistics  of  the  Discrete  Categories 
of  Refraction  and  Associated  Specifying  Parameters:  September 
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Statistics  of  the  Discrete  Categories 
of  Refraction  and  Associated  Specifying  Parameters:  October 
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Statistics  of  the  Discrete  Categories 
of  Refraction  and  Associated  Specifying  Parameters:  November 
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Statistics  of  the  Discrete  Categories 
of  Refraction  and  Associated  Specifying  Parameters:  December 
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Table  20 

Number  of  Profiles  Exhibiting  Different  Combinations 
of  Refractive  Conditions 
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The  one  profile  with  4 ducts  had  no  super-refractive  layers  and  four  sub-refractive  layers. 
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Parameter  S in  Adjacent  Latitude  Zones 


LATITUC: 

29  .CN 

. 5N 

NUN  6 

E? 

PROFILES 

29 

. CN 

-74 

• £K 

PE^CFNT 

N 

S fc 

s s 

SW 

DM 

DS 

TCT 

N 

SW 

c c 

OW 

DM 

ns 

N 

175 

0 • 

- V 

2 3 

21 

1J 

13 

222 

N 

54 

25 

7 

6 

0 

w 

sw 

32 

15 

12 

5 

4 

5 

77 

SW 

41 

24 

15 

6 

5 

e 

2 4 • 5 N 

c 2 

■a 

V 

6 

5 

4 

▲ 

3 

1 

1B 

2-*. 

2 c 

16 

4 4 

27 

e 

✓ 

C 

c 

3 7 • = N 

cw 

17 

c 

2 

3 

2 

J 

29 

37. 5N 

CW 

5 6 

17 

6 

10 

6 

L 

CM 

i 

2 

r 

M 

J 

J 

5 

ON 

60 

4 C 

r 

4 

0 

0 

c s 

C 

t 

0 

r 

r* 

J 

n 

J 

CS 

j* 

- 

p 

L 

c 

G 

r 

L /S T I TUCT 

24  ,=N 

-27 

. 2 N 

NUN  5 

F ° 

pRCFIlr  S 

24 

• 5N 

-.77 

. 5N 

PS 

PC  E K T 

N 

S W 

S S 

CW 

CM 

DS 

TCT 

N 

sw 

2 c 

OW 

O’T 

cs 

N 

265 

56 

25 

12 

2 

1 

26  7 

N 

74 

14 

6 

3 

3 

r 

sw 

56 

16 

4 

3 

1 

1 

51 

SW 

65 

1 5 

4 

3 

1 

1 

2 7 • 5 N 

ss 

17 

1C 

2 

4 

r, 

‘J 

24 

37.&N 

2 s 

5 3 

29 

6 

11 

C 

•» 

4i.  on 

CW 

1 1 

c 

2 

4 

J 

w 

22 

41.  ON 

CW 

5 G 

22 

c 

16 

f* 

■* 

CM 

r 

2 

1 

0 

3 

V# 

3 

CN 

0 

66 

2 2 

0 

p 

i c 

cs 

1 

C 

C 

c 

3 

0 

1 

CS.. 

o 

i 

- 

LATITUDE 

27  ,=N 

-4  1 

. CN 

NUNES'? 

PROFILES 

27 

. 5N 

-41 

. CW 

FEDCE 

NT 

N 

SW 

2 2 

CW 

DM 

DS 

TOT 

N 

SW 

2 c 

cw 

DM 

os 

N 

4 3 C 

77 

27 

i 3 

2 

4 

5 2 0 

N 

76 

14 

e 

2 

n 

4 

SW 

75 

1 j 

11 

7 

2 

66 

SW 

57 

16 

16 

10 

3 

r 

<*  1 . C N 

22 

11 

1 

2 

1 

G 

3 

16 

41.  uN 

2 2 

61 

16 

16 

5 

0 

C 

h 6 . 0 N 

CW 

15 

4 

p 

l 

0 

n 

2 C 

46. CN 

CW 

75 

20 

c 

C 

0 

. 

DM 

C 

r 

p 

F 

66 

U 

DN 

o 

o 

c 

0 

p 

c 

CS 

c 

0 

- 

CS 

> 

■j 

' 

J 

[ 

LATITUDE 

41  .CN 

-46 

. CN 

NUN  EE# 

PROFILES 

41 

. ? N 

-46 

. rN 

PERCENT 

N 

SW 

SS 

DM 

CS 

TOT 

N 

SW 

c 2 

cw 

CM 

DS 

N 

2^6 

■»6 

1G 

c 

T 

> 

341 

N 

6 3 

1j 

2 

? 

*. 

SW 

1 4 
. 4 

11 

1 

2 

1 

7 

u o 

SW 

64 

2 2 

6 

4 

2 

r 

46.  ON 

SS 

16 

4 

3 

1 

3 

F 

J 

26 

ho.CN 

SS 

69 

15 

11 

3 

F 

21. IN 

CW 

c 

2 

G 

c 

c 

'j 

14 

51.  CN 

CW 

*4 

21 

(. 

14 

r* 

o 

t* 

CM 

1 

C 

C 

r * 

0 

J 

1 

ONI 

r 1 

6 

n 

** 

C 3 

1 

c 

L 

F 

. u 

l. 

1 

r SI  'JO 

r 

c 

P 

r 

C 

- 1 35- 


I 


Table  27 


Comparison  of  the  Distributions  of 
Parameter  D in  Adjacent  Latitude  Zones 


LATITUCE 

34.5 

N 

N-  3 7 

SW 

N 

19 

7 

SW 

13 

16 

29  . . N S3 

13 

1 7 

3 *.  «t  N C <» 

19 

21 

CM 

3 

2 3 

CS 

1 

16 

LATITUCE 

37. 5 
N 

N-4  1 
S W 

N 

55 

3 C 

SW 

41 

28 

3 4 ,'j  N SS 

27 

22 

37. 5N  CW 

29 

37 

C- 

6 

c 

CS 

im 

4 

latituc: 

41.  2 
N 

N-4  6 
SW 

N 

122 

2 7 

SW 

o 8 

3 7 

3 7 .5  N S3 

4 1 

35 

41. CN  CW 

33 

t 6 

CM 

6 

6 

CS 

1 

2 

LATITUCE 

46.3 

N 

N-t  1 
SW 

N 

95 

33 

S W 

6 4 

26 

41. CN  SS 

31 

24 

46.  IN  CW 

24 

19 

5N  NUMBER  PROFILES 


£ c 

CW 

CP 

CS 

TCT 

2 

C 

£ 

K, 

25 

12 

14 

G 

1 

52 

i 3 

2 - 

e 

V 

s. 

73 

28 

4 7 

7 

1 

123 

14 

33 

15 

C 

5 8 

1C 

2 6 

1 3 

2 

7G 

.N 

MMcEF 

PROFILES 

SS 

CW 

CM 

CS 

TCT 

12 

15 

U 

« 

116 

19 

2 

4 

i.,9 

16 

2 o 

J 

89 

25 

41 

1 2 

4 

153 

3 

13 

6 

t 

42 

2 

2 

- 

*4 

:n 

NUMEEk 

PRCFILES 

SS 

CW 

CM 

CS 

TCT 

17 

1 

176 

26 

2 3 

2 

w 

1 1 o 

2 1 

13 

2 

"* 

116 

29 

15 

6 

z 

125 

4 

1 4 

•4 

24 

r, 

0 

« 

W 

1 

0 

7 

NU M c E R 

PROFILES 

S S 

CW 

CM 

r c 

TCT 

9 

< a 

2 

Z 

1 3 5 

12 

8 

3 

o 

no 

1 j 

6 

j 

1 

72 

12 

17 

1 

1 

74 

3i4.5K-37.tN  PERCENT 


N 

SW 

5 s 

CW 

CM 

CS 

N 

54 

2 Q 

c 

14 

c 

j 

J 

SW 

25 

7 1 

w W 

2 2 

19 

0 

4 

X 

2 5 . , N 

SS 

24 

22 

1 7 

27 

6 

w 

34. 5N 

CW 

15 

17 

2 2 

26 

5 

u 

DM 

3 

23 

1 4 

38 

15 

c 

OS 

4 

2 2 

14 

37 

38 

z 

37 

• 5N 

-4  1 

. CN 

PERCEN  T 

N 

SW 

c c 

CW 

CM 

rs 

N 

47 

25 

1 L 

16 

Q 

j 

SW 

27 

25 

1 7 

16 

3 

_ 

34. 5N 

SS 

22 

2 4 

1 7 

22 

4 

‘J 

37. 5K 

^W 

13 

22 

2 2 

2 5 

7 

2 

IM 

14 

21 

15 

7 2 

14 

3 

CS 

14 

28 

1 4 

21 

1 4 

7 

41 

• U N 

-46 

. r N 

PERCCNT 

N 

SW 

i S 

CW 

CM 

CS 

N 

69 

15 

c 

5 

'i 

* 

SW 

43 

22 

16 

14 

1 

c 

37. 5N 

SS 

24 

21 

1 7 

15 

1 

r 

4 1 . 0 N 

CW 

24 

24 

2 1 

14 

4 

l 

CM 

17 

17 

11 

41 

11 

CS 

14 

28 

r 

4 2 

14 

c 

46 

. j N 

-51 

• * N 

PERCENT 

N 

S W 

c c 

CW 

CM 

CS 

N 

61 

2 1 

5 

8 

1 

4 

SW 

63 

23 

1 r 

7 

9 

41. IN 

s s 

43 

3 3 

1 2 

8 

3 

1 

46  ,-N 

ow 

32 

25 

1 6 

22 

1 

1 

1 


Table  28 


Comparison  of  the  Distributions  of 
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Change  in  Distribution  of  Parameter  S 
in  12  Hours  from  00Z:  January-March 
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Change  in  Distribution  of  Parameter  S 
in  12  Hours  from  12Z:  January-March 
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Table  34 

Change  in  Distribution  of  Parameter  S 
in  12  Hours  from  12Z:  April-June 
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Table  35 

Change  in  Distribution  of  Parameter  S 
in  24  Hours  from  00Z:  April-June 
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Change  in  Distribution  of  Parameter  S 
in  12  Hours  from  12Z:  July-September 
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Change  in  Distribution  of  Parameter  S 
in  24  Hours  from  00Z:  July-September 
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Change  in  Distribution  of  Parameter  S 
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Change  in  Distribution  of  Parameter  S 
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Change  in  Distribution  of  Parameter  S 
in  24  Hours  from  00Z:  October-December 
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Table  44 

Change  in  Distribution  of  Parameter  S 
in  24  Hours  from  12Z:  October-December 
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Change  in  Distribution  of  Parameter  D 
in  12  Hours  from  00Z:  January-March 
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Change  in  Distribution  of  Parameter  D 
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Change  in  Distribution  of  Parameter  D 
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Table  55 

Change  in  Distribution  of  Parameter  D 
in  24  Hours  from  00Z:  July-September 
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Change  in  Distribution  of  Parameter  D 


in  24  Hours  from  12Z 
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Change  in  Distribution  of  Parameter  D 
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Change  in  Distribution  of  Parameter  D 
in  24  Hours  from  00Z:  October-December 
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CDC  3100  Refractive  Index  Format  (24  bits/word) 
Variable  Number  of  Records/Report 


1st  Record 

Word  1 

0001  octal 

2 

ship  number 

3 

date  (year*10000+month*  100+-day) 

4 

hour 

5 

latitude,  longitude  (lat*  100 (-long)  (whole  degrees) 

I Records  (I  is  i 

number  of  levels) 

Word  1 

0001  octal 

2 

Refractive  Index  N value  (J) 

3 

Refractive  Index  M value  (J) 

4 

height  meters  (J) 

K Records  (K  is 

the  maximum  of  the  number  of  traps,  number  of  super- 

refractive  layers , and  the  number  of  sub-refractive  layers . 

all  three  are  zero  there  are  no  records  of  this  type.) 

Word  1 

0001  octal 

2 

ship  number 

3 

latitude 

4 

longitude 

5 

date 

6 

hour 

7 

number  of  traps 

8 

number  of  super-refractive  layers 

9 

number  of  sub-refractive  layers 

10 

duct  bottom  (K) 

11 

trap  bottom  (K) 

12 

trap  top  (K) 

13 

trap  intensity  (K) 

14 

super-refractive  intensity  (K) 

15 

super-refractive  bottom  (K) 

16 

super-refractive  top  (K) 

17 

sub-refractive  bottom  (K) 

18 

sub-refractive  top  (K) 

19 

weather  1 

20 

weather  2 

21 

weather  3 

-171- 


Word 


CDC  6500  Refractive  Index  Format  (60  bits/word) 
One  Record/Report 


1 ship  number 

2 latitude  (lat  degrees*100+lat  minutes) 

3 longitude  (long  degrees*  100 Hong  minutes) 

4 date  (yr*10000  rmonth*I00  >-day) 

5 hour 

6 number  of  levels — NOL 

7 weather  1 

8 weather  2 

9 weather  3 

10  zero 

A height  (meters) 

B refractive  index  N value 

C refractive  index  M value 

Words  ABC  repeated  for  1=1,  NOL 

D number  of  traps  NT 

E number  of  super-refractive  layers  NS 

F number  of  sub-refractive  layers  NB 

G duct  bottom 

H trap  bottom 

O trap  top 

P trap  intensity 

Words  G H O P repeated  for  K = 1 , NT 

Q super-refractive  bottom 

R super-refractive  top 

S super-refractive  intensity 

Words  Q R S repeated  for  L = 1,  NS 

T sub-refractive  bottom 

U sub-refractive  top 

Words  T W repeated  for  M = 1 , NB 

V weather  1 

W weather  2 

X weather  3 


‘i 

4 


■■■  i 


fi 
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Standardized  Format  for  RI  Data  (60  bits/word) 


Each  report  is  one  record. 


1st  Word  - Bit 


2nd  Word  to  (1 
Bit 


0 

1 bit 

blank 

1 to  4 

4 bits 

year-  56 

5 to  8 

4 bits 

month 

9 to  13 

5 bits 

day 

14  to  15 

2 bits 

hour  (0,  00Z;  3,  12Z) 

16  to  23 

8 bits 

latitude  in  tenths  - 2 80 

24  to  31 

8 bits 

longitude  in  tenths  - 1200 

32  to  44 

13  bits 

ship  number 

45  to  47 

3 bits 

number  of  traps  (K) 

48  to  50 

3 bits 

number  of  super- refractive  (KS) 

51  to  53 

3 bits 

number  of  sub-refractive  (KN) 

54  to  59 

6 bits 

number  of  levels  (NOL) 

KNOL+D/2)  Word 

0 to  14 

15  bits 

N parameter  I level 

15  to  29 

15  bits 

height  meters  I level 

30  to  44 

15  bits 

N parameter  (1+1)  level 

45  to  59 

15  bits 

height  meters  (1+1)  level 

The  levels  are  followed  by  1 word  for  each  trap  then  1 word  for  each  super- 
refractive  layer  then  1 word  for  each  sub-refractive  layer  in  the  following 
formats : 


15  bits 

15  bits 

15  bits 

15  bits 

Traps 

intensity 

base  duct 

base  trap 

top  trap 

Super-refractive 

intensity 

blank 

base  layer 

top  layer 

Sub-refractive 

blank 

base  layer 

top  layer 

The  intensity  is  the  same  delta  N value  as  were  on  the  3100  and  6500  old 
formats . 

The  Rl  layer  data  is  followed  by  the  S parameter  and  then  the  D parameter, 
each  scaled  binary  thirty. 
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